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Principles  and  Applications  of  Radio 
Frequency  Impedance  Probes 


I.  HISTORICAL  INTRODUCTION 

In  1957  Yeung  and  Sayers1  reported  that  the  RF  impedance  between  two  elec¬ 
trodes  immersed  in  a  plasma  showed  a  minimum  at  a  certain  frequency  (Fig- 

**  2 

ure  1).  They  believed  it  was  the  plasma  frequency  u  .  In  1960  Takayama  et  al 
3  P 

and  Miyazaki  et  al  examined  a  Langmuir  probe  response  with  an  RF  modulation 

on  its  dc  voltage  and  found  a  sharp  increase  of  current  at  (Figure  2).  They, 

too,  believed  it  was  plpsma  frequency  w  . 

4  P 

In  1963  Levitskii  and  Shashurin  used  two  fine  wires  as  transmitter  and  re¬ 
ceiver  of  RF  signals  across  a  plasma  medium  and  discovered  that  the  resonance 


(Received  for  publication  22  July  1981) 

1.  Yeung,  T.  H.  Y.  and  Sayers,  J.  (1957)  An  RF  probe  technique  for  the  measure 

ment  of  plasma  electron  concentrations  in  the  presence  of  negative  ions, 
Proc.  Phys.  Soc.  70B:663-668. 

2.  Takayana,  K.  ,  Ikegami,  H.  and  Miyazaki,  S.  (1960)  Plasma  resonance  in  a 

radio -frequency  probe,  Phys.  Rev.  Lett.  5:238-240. 

3.  Miyazaki,  S. ,  Hirao,  K. ,  Aono,  Y. ,  Takayama,  K.,  Ikegami,  H.  and 

Ichiymma,  T.  (i960)  Resonance  probe  —  A  new  method  for  measuring 
electron  density  and  electron  temperature  in  the  ionosphere,  Rept. 
Ionosphere  Space  Res.  Japan  14:148-159. 

4.  Levitskii,  S.  M.  and  Shashurin,  J.P.  (1963)  Transmission  of  a  signal  between 

two  hi#i -frequency  probes  in  a  plasma,  Spy.  Phys.  Tech.  Phys.  54:319. 
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FREQUENCY  FREQUENCY 


Figure  1.  Yeung  and  Sayers'  Ex-  Figure  2.  RF  Modulation  on  the  dc 

periment:  Minimum  Impedance  at  Voltage  of  a  Langmuir  Probe 

Frequency  up 

5 

frequency  depended  on  sheath  thickness.  In  1964  Harp  showed  that  the  frequencies 
of  these  experiments  should  coincide,  depending  on  sheath  thickness,  and  be  slight¬ 
ly  below  the  plasma  frequency. 


5.  Harp,  R.S.  (1964)  An  analysis  of  the  behavior  of  the  resonance  probe,  Appl. 
Phys.  Lett.  4:186. 
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Consider  now  a  simple  parallel  plate  probe  system,  as  in  Figures  3  and  4. 
Theoretical  work  of  Harp  (1964)5  and  Crawford  and  Harp  (1965)^  suggested  that 
the  simple  plasma  probe  system  can  be  studied  as  an  equivalent  circuit  as  in 
Figure  4.  Ignoring  Lp  for  smallness,  we  have  a  series  of  Cg  -  Cp  -  Cg. 


where. 


C 


s 


A  A 
es  S  ~  S  ' 


(2) 


where  A  is  a  cross  sectional  area  of  probe,  and  eg(*  1)  is  the  dielectric  constant  of 
vacuum  sheath. 


C  = 
P 


A 

e  — 
P  P 


u  2  -  up2/(l  +  p/2S> 


2  2 
“  - 


so  that. 


Ceff(u)  "  Co 


U) 


-  CJ 


u2  -  up2/ +  p/2S) 


(3) 


where. 


C  = 
o 


e 


A 


s  d 


and 

d  -  p  +  2S 


Since  current  I(u)  is  given  by  I(w)  =  V(w)/Z(w), 

I(w)  =  ju  VM  C(u)  ,  (4) 

6.  Crawford,  F.  W.  and  Harp,  R.S.  (1965)  The  resonance  probe  —  A  tool  for 
ionospheric  and  space  research,  J.  Geophys.  Res.  70:587. 
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Figure  3.  Probe  System  With  Sheaths 


Figure  4.  Equivalent  Circuit 


For  the  system  considered. 


In  Eq.  (5)  the  current  I(u)  or  admittance  l/  |  Z(w)|  peaks  atu  =  uR  and  has  a 
zero  at  u  =  u^,  (Figure  5).  The  peak  is  finite  due  to  the  presence  of  resistance. 


WR  “p 

Figure  5.  Resonances  of  Current  in  Plasma  Condenser  System 

5 

Harp  (1964)  also  showed  that  for  a  spherical  RF  probe,  the  resonance  fre 
quency  Wp  is  given  by 


2 


1  +  R/(k  \D) 


(6) 


where  R  is  radius  of  probe,  AD  Debye  radius,  and  k  an  empirical  constant. 


3.  PLASMA  WITH  MAGNETIC  HELD 


In  the  presence  of  magnetic  field  B  the  system  becomes  anisotropic.  Charged 
particle  motion  parallel  to  B  is  unaffected  by  B  because  the  Lorentz  force 
e(V  X  B)  is  zero.  For  V  not  parallel  to  B,  helical  motion  results  (Figure  6). 

The  dielectric  function  e  becomes  a  matrix,  or  so-called  dielectric  tensor,  due 
to  the  anisotropy  of  space.  The  principal  modes,  parallel  to  or  perpendicular  to 
B,  are  relatively  simple 


E 


II 


where  u 

c 


—  is  called  cyclotron  frequency  (or  gyrofrequency). 


(7) 


B  FIELD  LINES 


Figure  6.  Helical  Particle  Trajectory 


It  is  known  that  the  frequency  of  plasma  oscillations  perpendicular  to  magnetic 
field  is 


UUHF 


(8) 
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which  is  the  frequency  when  e^(u)  =  0.  For  some  reason,  this  frequency  is  called 
the  upper  hybrid  frequency.  (The  lower  hybrid  frequency  is  for  ions.  )  The  plasma 
frequency  parallel  to  B  field  is  unchanged  because  p )  =  0. 

Consider  again  the  series  Cg  -  -  Cg  capacitance,  as  in  Figure  3.  With  B 

field  perpendicular  to  probe  electric  field  lines,  we  have  again 


_i_  =  J_  +  _L  +  J_ 

c  „  c  c  c 

eff  s  p  s 


where. 


C  *  - 
s  s 


2  2 
u  -  u 

c 


The  result  is: 


2  2 

-  u  ~  UUHF 

'o  2  ,  2  2S  2. 

u  -  (w„  +  — r  u  _  ) 
„  C  d  D 


where. 


Sheath  Hybrid  Resonance 


-o  '  ef 


d  =  p  +  2S  . 


In  the  limit  of  vanishing  magnetic  field,  we  find 


which  is  the  result  obtained  for  no-magnetic  field  case  in  Eq.  (3). 

4.  PHASE  AND  AMPLITUDE  OF  ADMITTANCE 

The  properties  of  the  complex  admittance,  derived  in  the  previous  section, 
are  now  investigated  in  more  detail.  The  admittance  is  (Figure  7) 

where, 

C  -  A 
Co  '  3 

d  =  p  +  2S  . 


2 

^SHR 


Figure  7.  Poles  and  Zeros  of  Admittance  of  RF 
Probe 


In  reality  there  is  always  energy  loss,  due  to  particle  collisions,  heating  of 
plasma  medium,  excitation  of  plasma  waves,  etc.  In  other  words,  there  is  a 
resistance  in  the  RF  probe.  This  fact  renders  the  infinity  of  the  amplitude  of 
admittance  to  be  finite. 

In  a  complex  plane  (Figure  8)  a  vector  A  represented  by  its  phase  and  ampli¬ 
tude,  defined  by 

Phase  6  =  tan'1  — i- |  (11) 

^Re  A  J 

and 

Amplitude  |A|  =  ^ /(lm  K)^  +  (Re  •  (12) 


Figure  8.  Amplitude  and  Phase  of  a 
Complex  Vector 


With  resistance  in  the  RF  probe,  the  behavior  of  the  admittance  of  the  probe 
is  shown  in  Figure  9.  It  is  remarkable  that  as  frequency  increases  beyond  sheath 
hybrid  resonance,  the  phase  changes  signature,  indicating  that  the  system  behaves 
now  as  an  inductance.  It  remains  behaving  inductively  until  the  frequency  increases 
beyond  WuHF*  ^  uPPer  hybrid  frequency.  Then  it  becomes  a  capacitance  again. 
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Figure  9.  Phase  and  Amplitude  of  Admittance  of  RF  Probe 


5.  TECHNIQUES  OF  USING  RF  PROBES  FOR  IONOSPHERIC  RESEARCH 

Some  of  the  noted  experiments  using  RF  probes  for  ionospheric  research  are 
now  reviewed.  Mackenzie  and  Sayers  used  a  parallel  plate  capacitor  on  Aeriel-3 
satellite  for  measuring  electron  density  at  various  altitudes.  Figure  10  shows 
a  simplified  diagram  of  their  instrument  circuits.  A  sawtoothed  voltage  sweep 
was  applied  on  the  capacitor.  When  the  capacitance  was  minimum  the  effective 
plasma  density  was  maximum  in  the  capacitor.  Thus,  the  sheath  effect  was  mini¬ 
mized  but  not  eliminated.  The  electron  density  was  deduced  from  the  plasma  fre¬ 
quency  obtained.  Figure  11  shows  the  electron  density  deduced  plotted  versus 
altitude.  The  general  behavior  of  the  graph  agrees  with  that  obtained  by  other 
methods  such  as  Langmuir  probe  method.  The  Langmuir  probe  method  is  known 
to  be  not  too  accurate.  Note  that  the  upleg  measurements  are  generally  higher 
than  the  downleg  measurements  for  some  unclear  reasons,  perhaps  outgassing 
from  the  rocket. 


7.  Mackenzie,  E.C.  and  Sayers,  J.  (1966)  The  radio  frequency  electron  density 
probe  for  rocket  investigation  of  the  ionosphere.  Planet.  Space  Sci.  14:731. 


/\y\  BIAS  VOLTAGE 


Figure  10.  RF  Probe  Circuit  Arrangement  With  dc 
Bias  to  Reduce  Sheath  Thickness  [  from  Mackenzie 
and  Sayers  (1966)^] 


Figure  11.  Electron-density  Profile  Obtained 
With  the  Birmingham  University  Probe  [  from 
Mackenzie  and  Sayers  (1966)'] 


In  later  experiments  short  cylindrical  antennas  or  spheres  were  often  used. 
Instead  of  applying  sawtoothed  voltage  to  minimize  sheath  effect,  the  upper  hybrid 
frequency  resonance  was  used.  The  sheath  hybrid  resonance  depends  on  the  shape 

of  the  capacitor  and  the  profile  of  sheath  density.  Therefore  it  is  not  usable. 

8 

Some  of  these  experiments  were  performed  by  Melzner  (1971)  on  Nike-Apache 
rockets,  Oya  and  Obayashi  (1966,®  1967, 10  1968**)  using  KAPPA  8,  KAPPA  9,  and 
L3H  rockets  and  Oya  and  Morioka  (1975)*^  on  TAIYO  Satellite.  All  of  these  experi¬ 
ments  used  the  upper  hybrid  frequency  resonance  technique.  An  accuracy  with 

±  1  percent  over  the  frequency  range  0.  7-17  MHz  in  0.  3  sec  has  been  claimed 

13 

[Ejiri  and  Obayashi  (1970)  ].  Examples  of  data  records  obtained  in  some  of 

these  experiments  are  shown  in  Figures  12  and  19.  Probable  ranges  of  electron 
and  magnetic  field  parameters  to  be  encountered  in  ionospheric  and  space  appli¬ 
cation  are  shown  in  Figures  14  and  15. 

Recently,  radio  frequency  impedance  probes  have  been  used  for  the  observa¬ 
tion  of  electron  beam-induced  effects  in  a  "mother-daughter"  rocket  configuration 
[  Maehlum  et  al  (1980)*“*]. 


6.  DISCISSION 


The  advantage  of  the  upper  hybrid  frequency  resonance  technique  is  its  inde¬ 
pendence  of  sheath  effect.  In  Eq.  (io)  the  sheath  effect  is  in  the  denominator 
while  the  numerator  is  unaffected.  It  is  the  numerator  that  governs  the  zero  of 
the  admittance.  It  is  true  that  the  zero  may  be  interfered  by  noise  level,  but  the 
phase  can  provide  an  accurate  zero  crossing  at  the  same  frequency 


8.  Melzner  (1971)  Plasma  Waves  in  Laboratory  and  Space.  J.  Thomas  (ed. )  2. 

9.  Oya,  H.  and  Obayashi,  T.  (1966)  Measurement  of  ionospheric  electron  density 

by  a  gyro-plasma  probe:  A  rocket  experiment  by  a  new  impedance  probe. 

Rep.  Iono.  Space  Res.  Japan  20:199. 

10.  Oya,  H.  and  Obayashi,  T.  (1967)  Rocket  measurement  of  the  ionospheric 

miasma  by  gyro-plasma  probe,  Rep.  lonos.  Space  Res.  Japan  21:1. 

11.  Oya,  H.  and  Obayashi,  T.  (1968)  Gyro-plasma  probe  measurement  of  the 

electron  density  in  the  transition  region  between  the  ionosphere  and  the 
magnetosphere,  Space  Res.  8:339-344. 

12.  Oya,  H.  and  Morioka,  A.  (1975)  Instrument  and  observations  of  gyro-plasma 

probe  installed  on  TAIYO  for  measurement  of  ionospheric  plasma  parameters 
and  low  energetic  particle  effects,  J.  Geomag.  Geoelectr.  27:331. 

13.  Ejiri,  M.  and  Obayashi,  T.  (1970)  Measurement  of  ionosphere  by  the  gyro- 

plasma  probe,  J.  Geomag.  Geoele.  1-12. 

14.  Maehlum,  B.N.,  Maseide,  K. ,  Aarsnes,  K.  ,  Egeland,  A.,  Grandal,  B.  , 

Holtet,  J. ,  Jacobsen,  T.  A.,  Maynard,  N.  C. ,  Soraas,  F. ,  Stadsnes,  J. , 
Thrane,  E.V. ,  Troim,  J.  (1980)  Polar  5  —  An  electron  accelerator  experi¬ 
ment  within  an  aurora.  1.  Instrumentation  and  Geophysical  Conditions, 
Planet.  Space  Sci.  28:259-278. 
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Figure  14.  Probable  Range  of  Electron  Density  and 
Temperature  to  be  Encountered  in  Ionospheric  and 
Space  Applications  l  from  Crawford  and  Harp  (1965)6} 
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MAGNETIC  FIELD  STRENGTH  GAUSS 


Figure  15.  Ranges  of  Application  of  the  Radio  Fre¬ 
quency  Probe,  wr  is  the  frequency  at  the  zero  of 
admittance,  and  r.  is  the  Larmer  radius  [from 
Crawford  and  Harp  (1965)6] 


19 


Furthermore,  it  can  be  shown  (Appendix)  that  the  numerator  of  Eq.  (10)  is  inde¬ 
pendent  of  the  profile  of  sheath  density,  as  long  as  the  sheath  distance  is  less 
than  d/2.  The  cyclotron  frequency  uc  can  be  measured  accurately  by  using  on¬ 
board  magnetometers.  Then  plasma  frequency  is  given  by  %/u^  -  to  uhk  • 

Not  only  the  frequency  of  the  zero  of  admittance  is  unaffected  by  the  density 
profile  of  the  sheath  around  a  probe,  but  also  it  is  insensitive  to  the  geometry  of 
the  sheath  and  probe. 

There  are  certain  limitations  on  the  capacitance  probe  methods.  The  capaci¬ 
tance  is  proportional  to  A,  the  area  of  capacitor  surface.  If  A  is  small,  the 
measured  value  of  capacitance  would  be  small.  This  fact  renders  the  use  of  capa¬ 
citance  probes  to  measure  the  plasma  density  in  a  small  area  practically  impossible. 

Another  limitation  is  the  time  needed  to  sweep  a  frequency  range.  Since  a 
rocket  or  space  vehicle  is  traveling  at  high  speed,  a  long  time  span  would  mean 
large  uncertainty  in  position. 
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Appendix  A 


Theorem:  The  Independence  of  the  Zero  of  Admittance 
on  the  Sheath  Electron  Density  Profile 


For  simplicity  and  clarity  of  discussion,  a  vacuum  sheath  has  been  used 
throughout  in  the  text.  It  has  been  shown  that  the  admittance  of  a  capacitor,  with 
a  vacuum  sheath,  has  a  zero  at  (or  Wyjjp  if  magnetic  field  is  present),  and  a 
pole  at  UgHR,  which  is  due  to  sheath  resonance.  While  the  pole  depends  on  sheath 
property,  the  zero  is  unaffected.  It  is  this  remarkable  property  that  renders  the 
use  of  the  radio  frequency  impedance  probe  an  accurate  method  for  determining 
plasma  density.  It  would  be  interesting  to  conjecture  that  this  remarkable  property 
remains  true  for  a  general  sheath  electron  density  profile;  it  will  be  shown  to  be 
true  in  this  appendix. 

We  generalize  the  probe  system  shown  in  Figure  3.  We  consider  a  sheath- 
plasma-sheath  in  series  (p  >  0,  s  >  0),  but  the  sheaths  are  no  longer  vacuum. 

Near  the  probe  walls,  electron  density  decreases  monotonically  with  the  distance 
from  the  wall.  The  effective  capacitance  Ce^  is  given  by 

i  r  .  i  .  i  <Ai> 

5  I  ds  p  /„»  +  p 

eff  J  ^sheath'3'  plasma 

where. 
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Cplasma 
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Csheath(s) 


^p2<s) 
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The  gyrofrequency  uc  is  independent  of  electron  density.  The  plasma  fre¬ 
quency  uip(s)  is  given  by 
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(s) 


4re2  N(s) 
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CA4) 


where  a  general  realistic  model  of  N(s)  is  shown  in  Figure  Al. 

N(s)  =  Nq  (1  -  exp  (-as))  ,  (A5) 

where  a  is  a  sheath  parameter  of  electron  density  profile. 


Figure  Al.  Sheath  Electron  Den¬ 
sity  Profile 


Consider  the  integral 


d 

1  f  ds  a  -  b  exp  (-as) 
o 


(A6) 


where  a  >  0,  b  >  0,  and  a  >  0.  Let  r  =  exp  (-as).  Integrating  over  r.  One 
obtains 


4 


2 


In  a  -  br 


2  2  2 

Substituting  a  =  Cq  and  b  *  CqNo  d  /m  (u  -  wc  )  into  Eq.  (A7)  and  using 


Eqs.  (Al),  (A2),  and  (A3),  one  obtains 


t  J  f  u2  -  uc2  2  u2  -  (Jc2  -  U  2  exp  (-od)  1 

C“7  *  2\  +  I  2  /  2  2T  I  +  2od)* 

eff  o  lu  -  (u  +  u  )  u  -  (ui  +  u  )  | 

c  p  c  p  ) 


where  the  factor  2  is  due  to  the  two  sheaths  in  a  parallel  plate  condenser  (Fig¬ 
ure  3). 

The  admittance  is  obtained  by  using  Eq.  (A8): 


2  2 
“  '  u  UHF 


a  c<  — 2771  2  : 

w  -  w  +  (w  -  UImE.) 


where 


2  2 
p  ,  u  '  uc 

F<“)  "  I - ~2 

Cu 


C  -  Up  exp(-od) 

£2 - 1 - 1  +  2ad  * 


From  Eq.  (A9)  it  is  readily  observed  that  the  zero  of  admittance  is  still  at  u  jjjfp 
while  the  pole  is  a  function  of  sheath  parameter  a. 

In  the  limiting  case  that  a  is  small  so  that 


Lim  F(o)  =  0 
o-*0 


The  pole  of  l/Z (w)  becomes  at  wc,  agreeing  with  the  limiting  result  in  Eq.  (10) 
with  no  sheath  (s  ■*  0  in  Eq.  (10)].  The  conjecture  is  proved. 


